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Abstract 

The weak boson fusion process for neutral Higgs boson production is investigated with particular 
attention to the accuracy with which the Higgs boson coupling to weak bosons can be determined 
at CERN Large Hadron Collider (LHC) energies in final states that contain a Higgs boson plus 
at least two jets. Using fully differential perturbative matrix elements for the weak boson fusion 
signal process and for the QCD background processes, we generate events in which a Higgs boson is 
produced along with two jets that carry large transverse momentum. The effectiveness of different 
prescriptions to enhance the signal to background ratio is studied, and the expected signal purities 
are calculated in each case. We find that a simple cut on the rapidity of one final-state jet works 
well. We determine that an accuracy of Sg/g ~ 10% on the effective coupling g may be possible 
after ~ 200 fb^^ of integrated luminosity is accumulated at the LHC. 
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I. INTRODUCTION 



Following the discovery of the neutral Higgs boson H at the CERN Large Hadron Col- 
lider (LHC), attention will focus on the measurement of its couplings to gauge bosons and 
fermions. A promising reaction from which to extract some of these couplings, particularly 
the HWW coupling, is the weak-boson fusion (WBF) process , where the Higgs 

boson H is produced via fusion of the weak bosons W and Z: WW, ZZ —>■ H, and is accom- 
panied in the final state by two jets that carry large transverse momentum pt- To extract 
the couplings reliably, a good understanding is required of the production processes and the 
background processes that lead to the same final state. Many strong interactions subpro- 
cesses also generate Higgs boson-plus-two-jet {H + 2 jet) final states. These background 
subprocesses can be computed with the techniques of perturbative quantum chromodynam- 
ics (QCD). They supply an irreducible background that may be reduced to some extent by 
judicious selections on the final state event topology. 

In the analysis presented here, we have in mind a situation in which the Higgs boson has 
been discovered and a sample of events exists containing a Higgs boson and two or more 
jets. This set of events will contain backgrounds of two types: real H + 2 jet events produced 
by QCD mechanisms other than WBF, and events which contain jets and particles that are 
present in typical Higgs boson decay modes, but without an explicit Higgs boson. Within 
the full event sample, we discuss the simulation of the real WBF signal and the irreducible 
QCD H + 2 jet background. We do not address the second type of contamination, such as 
events from the QCD Z + 2 jets process where the Z decay imitates a Higgs boson decay. 
Our concern is to estimate the expected signal purity, by which we mean the fraction of real 
Higgs boson events produced by weak boson fusion. 

The WBF H + 2 ]et signal region is characterized by jets that carry large transverse 
momentum and large rapidity. Because the jets carry large transverse momentum, it is 
necessary to use hard QCD matrix elements in order to represent the signal and the H + 2 iet 
background reliably. A parton shower approach, for example, would be unlikely to provide 
a correct estimate of the momentum distribution of the jets in the region of phase space 
of interest. Next-to-leading order (NLO) QCD^ corrections to the total WBF production 



cross section have been known for some time jl5 
calculated recently in a fully differential way [jj 



and the corresponding corrections were 
In this paper, we use an independent 
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calculation to verify the results of Ref. [16[ and to examine in more detail the effects of 
the WBF selection cuts on the NLO QCD corrections. We also use perturbative QCD 
expressions for the background H + 2 iet matrix elements. At present, the fully differential 
if + 2 jet background distributions are known only at leading order. In addition to our NLO 
study of the signal process, we nevertheless provide two estimates of the NLO enhancement 
of the QCD H + 2 jet background process, in order to better assess the viability of the 
WBF channel for measuring the coupling strength of the Higgs boson to vector bosons. Our 
calculations are fully differential at the partonic level. One limitation of the fact that we 
omit showering is that forward beam jets, which likely have low pr, are ignored. 

Since the WBF channel is most interesting for a Higgs boson in the mass range mn = 
115 — 200 GeV, we perform calculations with the two extremal values of this range. We 
compute differential cross sections in rapidity and transverse momentum at a pp collider with 
^/s = 14 TeV. In Sec. Ull we discuss the production processes that contribute to the WBF 
signal and backgrounds, and we describe our method for evaluating them. We gene rate 
momentum distributions using the general purpose Monte Carlo program MCFM '17]. In 
our case, all jets carry a minimum value of relatively large transverse momentum whose 
values we specify. 

We present numerical values of the differential cross sections in Sec. IIIII Various pre- 
scriptions are used in the literature to define the WBF sample, cuts that enhance the WBF 
fraction of the cross section by exploiting the special character of WBF events. Our inves- 
tigations lead us to propose a new, somewhat simpler definition in terms of a cut on the 
rapidity of one of the final state jets. In this section, we also define quantitatively what 
we mean by WBF signal purity P. We find that purities of 60% to 70% can be expected 
if a selection of pt > 40 GeV is made on the tagging jets and somewhat lower values if 
the cut is dropped to 20 GeV. We derive an expression for the expected uncertainty on the 
effective Higgs-boson-to-weak-boson coupling strength g in terms of P, the expected statis- 
tical accuracy of LHC experiments, and the uncertainties on the signal and the background 
processes. We estimate that it should be possible to achieve an accuracy of Sg/g ~ 10% 
after ~ 200 fb^^ of integrated luminosity is accumulated at the LHC. Somewhat smaller 
values of Sg/g are obtained in another recent investigation of anticipated uncertainties in 
the couplings and we explain the source of the difference. 

In Sec. IIV| we compare the effects on both event rates and signal purity of our proposed 
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method for defining WBF events with two other popular methods found in the literature: 
a selection on the difference in rapidities between two tagging jets in the final state, and 
a selection on the invariant mass of a pair of tagging jets. The alternative prescriptions 
yield some increase in signal purity with respect to our definition, but the gain is sensitive 
to the cut in transverse momentum used to specify the trigger jets, and it is accompanied 
by loss of event rate. For values of the jet cut pr > 40 GeV, our prescription appears to 
work about as well as the other methods. Relatively high luminosity will be needed for a 
precise determination of the uncertainty Sg/g. Our simpler definition of the WBF sample 
in terms of a selection on the rapidity of only one jet offers advantages in a high luminosity 
environment where a large value of the transverse momentum cut is appropriate and multiple 
events per crossing may be an issue. 

We provide two methods for estimating the size of next-to-leading order corrections to the 
H + 2 jet background in Sec. |3 One of these relies on similarity with the Z + 2 jet process for 
which fully differential NLO results are known. The other method is an extrapolation from 
the known next-to-next-to-leading order (NNLO) results for the fully inclusive Higgs boson 
production process. The substantially different estimates for the NLO enhancement provided 
by these two methods show the level of uncertainty of the LO background calculation. A 
fully differential NLO calculation of the if + 2 jet background applicable in the region of 
interest for WBF investigations is needed in order to improve our computations of signal 
purity and of the expected uncertainty in 6g/g. A summary of our conclusions may be found 
in Sec. ED 

II. PRODUCTION CROSS SECTIONS 

Examples of the WBF diagrams that must be calculated are shown in Fig. ^ The basic 
leading order process is shown in (a), where the exchanged bosons may be either W^s or Z's, 
and one or both quark lines may be reversed, yielding qq, qq and qq initial states. The virtual 
NLO corrections are obtained by adding a gluon loop to either qqV vertex, as illustrated in 
(b). The remaining real NLO corrections are shown in (c), where either an additional gluon is 
radiated in the final state or a gluon from the proton splits into a qq pair. Calculation of the 
necessary loop diagrams is straightforward, providing a couple of simplifying assumptions 
are made. First, we ignore contributions of the form qq' —>■ V* —>■ VH, where V = W, Z. 
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(a) 



(b) 




(c) 



FIG. 1: Representative diagrams for the production of a Higgs boson via weak boson fusion: (a) at 
fowest order; (b),(c) at NLO. Further diagrams can be obtained by crossing incoming and outgoing 
lines in all cases. All of the virtual corrections are of the vertex correction form, as shown in 
(b). There are two types of real corrections depicted in (c). The first set corresponds to the 
emission of a gluon in all possible positions on the quark lines (left-hand diagram) and the second 
set corresponds to the crossing where a gluon is present in the initial state (right-hand side). 

Second, we neglect any interference effects from identical flavor quarks in the final state. We 
checked that both of these approximations have little effect on the calculated cross sections 
at leading order, particularly in the region of phase space that we consider. 

The NLO calculation is embedded in the general purpose Monte Carlo program 
MCFM which uses the dipole subtraction method 3]. We use the default set of 
parameters in this program, in which a = 1/128.89, Mw = 80.419 GeV, Mz = 91.188 GeV, 
and sin^ 6'^ = 0.2285. For the parton distribution functions, we use CTEQ6L1 for lowest 
order and CTEQ6M at NLO j20|. In these sets of parton densities, a^°(M^) = 0.130 and 
af^'^lMz) = 0.118. In this paper, we choose the reference value fi = itlh for the renormal- 
ization and factorization scales. In Sec. IIII Cl we investigate the uncertainty of the signal and 
of the background associated with variation of the scale over the interval 2mH > /i > mul'^- 
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A. Generic cuts 



The hallmark of WBF events is a Higgs boson accompanied by two "tagging" jets having 
large pt and large rapidity. In real events and in computations at NLO, there are generally 
more than two jets, and the goal is to pick out a clean signal. To simplify our study and 
to demonstrate the robust character of the WBF process, it is desirable to make as few 
selections (cuts) as possible on the events. We begin with a basic set of cuts, exactly as 



in Ref. 



. Jets obtained from the Monte Carlo runs are clustered according to the kr 



algorithm with pj^ > 20 GeV, jet pseudo-rapidity \ri^°^\ < 4.5, and jet separation AR 
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A.r]jj + A.(f)jj > 0.8, where A0jj is the difference in the azimuthal angles of the two jets 
in the transverse plane. The two jets with the highest pt are chosen as the tagging jets 
and ordered according to their pseudo-rapidities, rjj-^ < rjj^. In order to approximate the 
acceptance for the Higgs boson decay products and to complete the specification of our 
minimal set of cuts, we imagine the decay of a Higgs boson to two charged particles, denoted 
as "leptons". We require that these leptons satisfy the cuts: 

p'^""' > 20 GeV, |r/i^P*| < 2.5, Ai?,, > 0.6, rjj, < r/iept < 

The Higgs boson decay products are therefore located in pseudo-rapidity between the two 
high-pT jets. Although we enforce these cuts on potential Higgs boson decay products, the 
cross sections that we present do not include any branching ratio for this decay. We include 
branching ratios and efficiencies when we discuss the determination of the coupling strength 
in Sec. IhTI 

Throughout this paper we refer to the QCD production of a Higgs boson in association 
with jets as the "background" to our WBF signal events. This cross section is implemented 
in MCFM at leading order based on the matrix elements of Ref. [2^. A selection of the 
contributing diagrams is shown in Fig.|21 There are contributions from qq, qg, and gg initial 
state subprocesses, but the Higgs boson is always produced from an effective gg ^ H vertex. 
The effective coupling of the Higgs boson to two gluons is included in the limit of heavy 
top-quark mass m^, with a coupling strength ^^/(Svrt'), where v is the Higgs boson vacuum 
expectation value, and as is evaluated at the scale rriH with the 1-loop expression for the 
evolution (0.125 for mn = 115 GeV and 0.116 for mn = 200 GeV). The effective coupling 
approximation should be valid since we limit ourselves to Higgs boson masses rriH < 2mt 
and Higgs boson transverse momenta pf^ < rrit 



221. 



FIG. 2: Representative diagrams for the production of a Higgs boson and two jets at lowest order, 
calculated in the heavy top-quark limit of the Hgg effective coupling. 



In Sec. 13 we estimate the NLO corrections to the lowest order result for the H + 2 jets 
back grou nd process by comparison with the similar Z + 2 jets process, calculated at NLO in 
Ref. [iTI, l23|. and by extrapolation from NNLO calculations of the fully inclusive process 



pp HX f2J-j30(, and NLO calculations oi pp ^ H + I ]ei + X 



III. RESULTS 



The cuts mentioned in the previous section are a generic set of cuts. They do not exploit 
the kinematic structure of WBF events, where the jets tend to be produced very forward 
in pseudo-rapidity. In this subsection, we present first the cross sections for the signal and 
background processes after application of the generic cuts. Without further cuts, the WBF 
events would be lost in the QCD continuum background. We then apply one additional 
constraint which defines our WBF sample, and we show results for kinematic distributions, 
event rates, and signal purities. 



A. Basic cuts 



We examine the effects of the generic cuts in terms of their effects on the WBF signal 
and the if + 2 jet background. These cross sections - without any further cuts - are shown 
in Table H] as a function of the minimum jet pt- The WBF signal process is calculated at 
NLO and, at this point, the if + 2 jet background at LO. We remark that the effects of 
the NLO corrections on the WBF process are rather small in this region, corresponding to 
i^-factors between 0.95 and 1.1, depending on the Higgs boson mass and px cut. 

The values in the table show that (without consideration yet of NLO effects in the back- 
ground) the rates for the signal and background are comparable for p™™ = 20 GeV, and the 
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Pt cut [GeV] 


20 


40 


80 


Signal {rriH = 115) 


1866 


1081 


239 


Bkg 


2173 


743 


200 


Signal {mn = 200) 


1189 


709 


166 


Bkg 


958 


340 


96 



TABLE I: Cross sections in fb for the WBF signal (calculated at NLO) and + 2 jet back- 
ground(LO), as a function of the minimum jet px- Only the minimal set of cuts of Sec. ITll is 
applied. 

signal-to-background ratio improves as the pt cut is increased. The WBF signal lies above 
the if + 2 jet cross section if p™" > 40 GeV. 

B. WBF cuts 

In an attempt to exploit the WBF event structure, a popular cut invokes a separation in 
pseudo-rapidity between the two tagging jets, for instance \rij-^ — rjj^ \ > 4. In this paper, we 
define a slightly different and simpler cut, motivated by our examination of the distributions 
of the absolute jet pseudo-rapidities shown in Fig. El In these figures, each tagging jet enters 
with weight one-half and cross sections have been converted to event rates with an integrated 
luminosity of 1 fb~^. The area under each curve is equal to the total number of events in 
that channel. 

The plots in Fig. IHl show that the shape of the distribution depends little on either the 
Higgs boson mass or the jet pt-, but - as expected - is very different in the signal process, 
compared to the if + 2 jet background. In each case, the WBF events peak at values of 
|?7| ~ 3, although there is a slight movement to lower values of |?7| as the pr cut is increased. 
The width of the peak also tends to decrease, but the full width at half-maximum is fairly 
constant at approximately 3 units of rapidity. In contrast, the rate of background events 
falls off fairly sharply beyond |?7| ~ 2. 

Motivated by the comparison of rapidity spectra in Fig. El and erring on the side of 
simplicity, we choose a uniform cut that ensures at least one jet lies within the peak. Namely, 

??peak - ^?width/2 < |?7il < r?peak + '7width/2, (1) 
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FIG. 3: Dependence of the tagging jet pseudo-rapidities on the minimum jet pT used, for the two 
cases TJiH = 115 GeV (left) and mn = 200 GeV (right). Each of the two tagging jets in the event 
is entered in these plots, with weight one-half, and the rates assume an integrated luminosity of 
1 fb~^. The signal (solid) is calculated at NLO and the background (dashed) at LO. We show 
results for three different selections of the minimum jet transverse momentum, px > 20, 40, and 
80 GeV. 
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Pt cut [GeV] 


20 


40 


80 


Signal {niH = 115) 


1374 


789 


166 


Bkg 


1196 


382 


92 


Purity 


0.53 


0.67 


0.64 


Signal {mn = 200) 


928 


545 


121 


Bkg 


534 


179 


46 


Purity 


0.63 


0.75 


0.72 



TABLE II: Event rates for the Hjj WBF signal(NLO) and Hjj background(LO), including our 
WBF requirement that at least one jet carry large \ri\, as defined by Eq. (pQ)- We assume 1 fb~^ 
of integrated luminosity. Purity is defined as P = S/{S + B), where S stands for the number of 
signal events and B for the number of background events. 

for j = ji or j = j2, where 77peak=3 and ?7width=2.8. 

Equation (OP, along with the generic cuts specified above, constitutes our definition of 
weak boson fusion cuts. The effects of the pseudo-rapidity restriction on the jets are shown 
in Table ITTl where we have assumed 1 fb~^ of integrated luminosity. The rates in this table 
should be contrasted with those in Table lU In Table |nj we include values for the signal 
purity, defined as P = S/{S + B), where S stands for the number of signal events and B 
for the number of background events. The number of events as a function of the minimum 
jet Pt is also plotted in Fig. |3] Comparing the tables, one can see that the signal rate is 
diminished only slightly, by about 20-30%. On the other hand, the background is shrunk 
considerably, by about a factor of two. A px cut of 20 GeV is barely sufficient to distinguish 
the WBF signal above the QCD LO Hjj background for ttlh = 115 GeV. However, the 
signal 5* to background B ratio improves to about 2 for p^^ > 40 GeV. At rriH = 200 GeV, 
the situation is better, with S/B of about 1.7 when the pt cut is 20 GeV, and rising to 
~ 3 for p™^ > 40 GeV. A px cut of 40 GeV yields a prominent effect across the range of 
interesting masses, ran = 115-200 GeV. 

It is instructive to examine the origin of the different rapidity spectra for the signal and 
the background. Since if + 2 jet events are generated in both cases with identical cuts on the 
transverse momenta of the jets, the different rapidity spectra must originate from dynamics. 
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FIG. 4: Numbers of events for the WBF signal and the QCD background as a function of the 
minimum jet px, for an integrated luminosity of 1 fb^^. No branching ratios for the Higgs boson 
decay have been applied. The pseudo-rapidity restriction on one of the jets has been enforced, as 
in Table Ull The solid line is the NLO signal and the dashed is the LO background. 

Comparing the LO production diagrams in Figs.C]and|21 we note that gg and qg initial states 
contribute to the QCD background but not to the WBF signal. The gluon parton density is 
notably softer than the quark parton density, suggesting a plausible reason for the differences 
in the rapidity spectra of the final state jets in the two cases. This reasoning is supported 
by the results shown in Fig. El The shape of the background rapidity spectrum from the qq, 
qq, and qq contributions is very similar to that of the signal, albeit with a slight shift of the 
peak to smaller \ri\. The very different rapidity spectra of the signal and the background 
evident in Fig. El results therefore primarily from the gg and qg initial state contributions. 
The results shown in Fig. El imply that there is a basic upper limit to the purity one can 
achieve for the WBF event sample, regardless of which prescription one adopts to define the 
WBF sample. The qq, qq, plus qq component of the QCD background process generates 
a final state event topology essentially identical to the WBF signal process. Values of the 
purity are listed Table ITTTl there is not much variation with tjih oi the value of the cut in 
Pt- Our results suggest that purity is bounded from above by at most P < 0.95 at LHC 
energies. 



C. Scale dependence study 



To examine further the effects of NLO corrections, we consider variation of the renormal- 
ization and factorization scale. A range of values < < 2m h is used conventionally 
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FIG. 5: Dependence of the tagging jet pseudo-rapidities for jet pT > 40 GeV, for the two cases 
niH = 115 GeV (left) and mn = 200 GeV (right). For the background, we show the full result 
with all contributions included and, for comparison of shapes, the background obtained if only the 
qq, qq, and qq initial state contributions are used. The magnitude of the separate component is 
multiplied by 20. 



PT cut [GeV] 


20 


40 


80 


Signal {rriH = 115) 


1374 


789 


166 


Bkg {qq,qq,qq) 


98 


45 


15 


Purity {qqmm) 


0.93 


0.95 


0.92 


Signal {mH = 200) 


928 


545 


121 


Bkg 


47 


23 


8 


Purity 


0.95 


0.96 


0.94 



TABLE III: Event rates for the Hjj WBF signal(NLO) and for the part of the Hjj back- 
ground(LO) that arises from the qq^ qq, and qq initial-state terms. 

to estimate the theoretical uncertainty at a given order in perturbation theory. As a rep- 
resentative case, we show results for a minimum jet pt of 40 GeV and both Higgs boson 
masses. In Fig. we show the tagging jet pseudo-rapidity distributions for the signal and 
background for a range of values of the common renormalization and factorization scale 
/i. The signal process shows very little variation with /i, a shift of less than ±2% when 
\rij\ ~ 2 in the WBF signal region. In contrast, the 17-1-2 jet background at LO is enhanced 
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FIG. 6: Tagging jet pseudo-rapidity distributions in Hjj events, calculated with a range of values 
of the renormalization and factorization scale, /x. The signal (solid) is calculated at NLO and the 
background (dashed) at LO. 

by approximately 70% when the scale choice fi = rnHfi is made, and reduced by 40% for 
fi = 2mH- A fully differential NLO calculation of the if + 2 jet background process is 
required to reduce the large uncertainty associated with fi variation apparent in Fig. IHl 

D. Uncertainty on the Couplings 

The signal and background events both include a real Higgs boson along with two jets. 
We may define a signal "purity" as the ratio S/ (S+B), where S denotes the number of signal 
events and B the number of background events. The purity as defined here does not improve 
with greater luminosity nor does it depend on the Higgs boson decay mode considered. Of 
interest to us is the effect of signal purity on the accuracy of the determination from data 
of the Higgs boson couplings gww and gzz to the WW and ZZ channels. The WBF cross 
section is proportional to a combination of g^^/ and g'^^, and their relative contribution 
changes somewhat with the value of the cut on pt- In this paper, we discuss only an effective 
coupling strength g. We remark also that in our discussion of the expected accuracy on g, 
we limit ourselves to uncertainties at the level of production of the Higgs boson. We set 
aside uncertainties associated with the fact that the Higgs boson is observed only in specific 
final states and that all the final states cannot be observed above backgrounds. 

To derive the uncertainty Sg/g on the coupling, we begin with the observed number of 
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events N = S + B. We define the ratio r = fi'observed/S'predictcd- Then, under the assumption 
that any deviation in the expected total number of events arises from the effective coupling, 
we obtain r = {N — B)/S. Taking the total derivative, we obtain an expression for the 
uncertainty in r. 

5r/r = ^[{5S/Sf + {{5Nf + {5Bf)/{N - Bf], (2) 
and, correspondingly, 

5g/g = l/2^[i6S/Sy + {{6Ny + {6By)/{N - B)^]. (3) 
With purity P = S / {S + B) , we derive 

6g/g = 1/2^[{5S/Sf + {l/Pf{5N/NY + ((1 - P)/ PfidB/ Bf]. (4) 

In the absence of any uncertainty in knowledge of the signal and background, Eq. Q 
demonstrates the obvious fact that the best one can achieve is 6g/g = 0.5 6N/N for P = 1. 
The factor (1/-P) that multiplies 6N/N in Eq. shows that reduction in the purity effec- 
tively reduces the statistical power of the data. Similarly, the factor (1— P) / P that multiplies 
6B/B shows that greater purity diminishes the role of uncertainty in our knowledge of the 
background. Given that purity decreases as the background increases, we see that the size 
of the background in the WBF region is the problem to contend with; the uncertainty on 
the background is of less importance. To represent the background reliably in a region of 
phase space in which tagging jets carry large transverse momentum, it is clearly important 
to use partonic hard matrix elements that can simulate this jet activity. 

To obtain numerical values for the overall uncertainty in g, we must specify the uncer- 
tainties expected in our knowledge of the signal S and the background B, along with the 
statistical uncertainty 6N/N. We address each of these contributions. 

The NLO QCD effects on the WBF signal are modest as are uncertainties associated 
with parton densities and variation of the renormalization/factorization scale. We may take 
6S/S = 5%, based on the calculated NLO /i dependence of ±2% and PDF uncertainty 
of ~ 3%, both obtained in the WBF region of phase-space. Next-to-leading order QCD 
contributions to the background are discussed in Sec. El The size of these contributions, 
scale dependence, and parton density variation could make the LO background estimate 
uncertain at the 60% level. We adopt a perhaps optimistic value of 6B/B = 30%. This 
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FIG. 7: The predicted uncertainty Sg/g in the couphng of the Higgs boson to a pair of W bosons 
is shown as a function of signal purity P = S/{S + B) for expected statistical accuracies 6N/N of 
10% and 2% The uncertainties in knowledge of the signal S and background B are assumed to be 
5% and 30% respectively. 

choice presupposes that the 20% fi variation and 5% PDF uncertainty of the fully inclusive 
NLO cross section for Higgs boson production may also apply to the NLO calculation of 
H + 2 ]et production in the WBF region of phase-space, once this calculation is done. 

Based on a study of conventional backgrounds 3^, a minimum of roughly 10 fb~^ of 
integrated luminosity is needed to discover the Higgs boson in the WBF process. This 
figure would be achieved after one year of LHC operation at a luminosity of 10^^cm~^s~^. 
Using the numbers in our Table |Hl we expect a WBF sample (5* + 5) of ~ 12000 events 
for niH = 115 GeV and pTcut = 40 GeV, and ~ 7000 events for rriH = 200 GeV and 
PTcut = 40 GeV. To translate these event rates into statements about statistical significance, 
we must specify a Higgs boson decay mode and approximate tagging efficiencies for the 
decay products. For rriH = 115 GeV, we choose the decay H t^t^ , with one r decaying 
leptonically and the other hadronically These choices yield a branching ratio of 

Br{H tt) X Br{T leptons) x 5r(r hadrons) = 0.073 x 0.7 x 0.65 = 0.033. 
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For the efficiency for tagging hadronic r decays, we take the figure 0.26 from Ref. j33| as an 
optimistic upper bound. ^ The true value for this efficiency will be known only after analysis 
of data from LHC experiments. While our choice of tagging efficiency may seem large, it 
is relatively easy to scale our final results if a different value is preferred. The combination 
of branching fraction and efficiency results in a reduction in the number events by a factor 
e = 0.033 X 0.26 ^ 0.01. For rriH = 200 GeV, the decay H W~^W~ is prominent, and we 
select the case in which both W^s decay leptonically |ll|.^ We obtain 

e = Br{H WW) x Br{W leptons)^ = 0.74 x 0.22^ = 0.036. 

Using these numbers, we compute expected statistical uncertainties of 6N/N ~ 10% and 
~ 6% at ruH = 115 and 200 GeV, respectively. With statistical accuracy 6N/N of 10%, 
SS/S = 5% and 5B/B = 30%, we obtain 5g/g - 10% for purity P = 0.7 when uih = 
115 GeV, and 5g/g ~ 8% when mn = 200 GeV. 

After 5 years of LHC operation, we can anticipate an integrated luminosity of 200 fb~^ 
will have been accumulated. This increase allows us to reduce our estimates of 6N/N to 
~ 2% and ~ 1.5% at tuh = 115 and 200 GeV, respectively, and 6g/g ~ 7% for P = 0.7. 

We remark that the uncertainties in the signal 5* and background B dominate the uncer- 
tainty in g. If P = 0.7 and 6N/N = 2%, then the uncertainties 6S/S and 6B/B would have 
to be reduced to 3% and 6%, respectively, before the statistical uncertainty would control 
the answer. Even if P = 1, 6g/g is controlled by 6S/S until 6S/S < 6N/N. 

In Fig. [71 we show numerical predictions for the uncertainty as a function of purity, for two 
choices of the statistical uncertainty. Signal purities of 0.65 or greater permit determinations 
of 6g/g of 10% or better after 200 fb~^ have been accumulated. As shown in Table ITU 
P > 0.65 is obtained for prcut > 40 GeV at mn = 115 GeV, and prcut > 20 GeV at 
rriH = 200 GeV. The curves indicate to us that there is not much to gain from purities 
greater than 70%. 

^ We recall that the generic acceptance cuts defined in Sec. m for the Higgs boson decay products are 
included in our event rates. 

^ We acknowledge that H W'^W^ with leptonic decay of both W^s is not a perfect match to our earlier 
specification of two body decay of the Higgs boson to "leptons" . Because the W is fairly massive, the 
rapidity distribution of the decay leptons may extend beyond |?7| < 2.5, and a further acceptance correction 
may have to be applied. Such a study is best addressed by an experimental simulation of the entire decay 
chain. 
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FIG. 8: The predicted uncertainty Sg/g in the couphng of the Higgs boson to a pair of W bosons 
is shown as a function of signal purity P = S/{S + B) for expected statistical accuracies 6N/N of 
10% and 2% The uncertainties in knowledge of the signal S and background B are assumed to be 
4% and 20% respectively. 

Somewhat smaller values of SS/S and ^BjB are chosen in another recent investigation 
of anticipated uncertainties in the couplings [1^. These values are 6S/S = 4% and 6B/B = 
20%. Although the scope of that study is quite different from ours, we may compare our 
estimates with theirs. In Fig. |H1 we show the uncertainty as a function of purity for these 
new estimates of 6S/S and 6B/B. For P = 0.7, we now find 6g/g ^ 9% and ~ 5% for 
the low- and high- luminosity data samples. This new lower value of 6g/g is similar to that 







obtained in Ref. [1^ at comparable luminosity 



One must bear in mind that the uncertainty discussed in Ref. [isj is the uncertainty on and therefore 
a factor of 2 greater. 
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IV. ALTERNATIVE DEFINITIONS OF THE WBF SAMPLE 



In the previous section, we define the WBF sample by a simple selection on the rapidity of 
one jet in an event in which there is a Higgs boson and two jets each carrying pjn greater than 
a specified minimum value. Other definitions have been used in the literature, and we wish 
to compare our signal rates and purities with those obtained if we use these alternatives. 
We examine the traditional cut on rapidity separation between the two trigger jets and a 
cut on the invariant mass of the pair of trigger jets. 

A. Rapidity separation cut 

In Fig. El we illustrate the expected event rates as a function of the difference in rapidities 
between the forward and backward tagging jets. There is a clear separation in the locations 
of the peaks of the WBF signal and the background, not unlike that seen in our Fig. 01 
Distributions such as these may motivate the choice of a cut on rapidity separation, \riji — 
T]j2\ > 4, as in Refs. 

Q H, and Q- 

Signal and background rates for WBF events selected 
in this fashion are shown in Table HVl Comparison of Tables IITI and II VI shows that the signal 
rate is diminished somewhat and that the purity is greater when the rapidity separation 
selection is made. As shown in Figs. |7| and |H1 a gain in purity reduces the uncertainty in 
6g/g. The quantitative shift from P = 0.67 to P = 0.78 at M = 115 GeV and pt > 40 GeV 
is an improvement of only 3% in 6g/g, and this reduction is offset somewhat by the loss in 
statistical accuracy. 

At lowest order in perturbation theory, one might expect naively that our simple rapidity 
selection and the rapidity separation cut are close to identical since there are only two jets 
in the event, tending to be located in opposite hemispheres. However, the finite rapidity 
carried by the Higgs boson introduces differences. Our preference for the simple rapidity 
cut is based on a few considerations. In data (and at yet-higher orders in perturbation 
theory), there will be many jets, and the simple specification of events that satisfy Eq. 
wiU be easier to implement. Second, in a high luminosity environment with more than one 
event per beam crossing, a selection on only one jet (in addition to the Higgs boson) reduces 
the chance that jets from different events are used. Finally, in our study of NLO event 
topologies with three jets in the final state, we find that a gluon jet, rather than a quark jet 
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FIG. 9: The difference between the two tagging jet pseudo-rapidities for a minimum jet px of 
20 GeV, for the two cases rriH = 115 GeV (left) and mn = 200 GeV (right). The rates assume 
an integrated luminosity of 1 fb~^. The signal (solid) is calculated at NLO and the background 
(dashed) at LO. 



PT cut [GeV] 


20 


40 


80 


Signal {rriH = 115) 


1297 


718 


137 


Bkg 


758 


207 


38 


Purity 


0.63 


0.78 


0.78 


Signal {mn = 200) 


911 


521 


106 


Bkg 


349 


102 


20 


Purity 


0.72 


0.84 


0.84 



TABLE IV: Event rates for the Hjj WBF signal(NLO) and Hjj background(LO), without our 
WBF definition and instead with the rapidity separation cut |r/j^ — r]j^ \ > 4. We assume 1 fb^i of 
luminosity. The purity and significance are as defined before. 

is sometimes one of the two jets with largest pt- For example, with a jet cut of 20 GeV and 
TJiH = 115 GeV, a gluon is a tagging jet about 25% of the time when we use our definition 
of WBF events. 
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B. Invariant mass cut 



As an alternative to the rapidity separation cut, one might consider a cut on the invariant 
mass of the two trigger jets. In Fig. and in Table |3 we display the the effects of the 
mass cut Mjj > 800 GeV. Comparison of Figs. [TUl and El shows a decided improvement in 
the signal to background ratio, an effect that is borne out in the purity numbers shown in 
Tables fVl and ITTl However, the significant gain in purity is true only for the smaller values of 
the pt cut and is accompanied by a substantial loss of signal rate. Since the smallest value 
of the Pt could be employed only with low- luminosity data samples, it is not evident that 
the price in loss of signal rate is affordable. The combination of the mass cut and our simple 
forward jet cut improves purity only slightly and reduces the signal rate further. 




FIG. 10: The tagging jet pseudo-rapidity distribution for a minimum jet px of 20 GeV, for the two 
cases rriH = 115 GeV (left) and mn = 200 GeV (right) - with a minimum dijet invariant mass, 
Mjj > 800 GeV. Each of the two tagging jets in the event is entered in these plots, with weight 
one-half, and the rates assume an integrated luminosity of 1 fb^^. The signal (solid) is calculated 
at NLO and the background (dashed) at LO. 

Using the invariant mass cut to define the WBF sample, we note that the number of 
events at rriH = 115 GeV with a cut on pt of 20 GeV is very similar to what we obtain 
with our WBF prescription but with a cut on pj- = 40 GeV. The purities are also nearly 
the same. This comparison would seem to favor our simple prescription: a larger value 
of the cut on pt is more appropriate in a high rate environment and more effective at 
reducing backgrounds not considered here. We conclude this section with the remark that 
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Pt cut [GeV] 


20 


40 


80 


Signal {rriH = 115) 


808 


561 


158 


Bkg 


304 


183 


82 


Purity 


0.73 


0.75 


0.66 


Signal {niH = 200) 


617 


428 


121 


Bkg 


157 


95 


43 


Purity 


0.80 


0.82 


0.74 



TABLE V: Event rates for the Hjj WBF signal(NLO) and Hjj background(LO), without our 
usual WBF definition and instead with Mjj > 800 GeV. We assume 1 fb^^ of luminosity. 



the alternative prescriptions of the WBF sample in terms of either a rapidity separation cut 
or an invariant mass cut yield some increase in the signal purity with respect to our simple 
cut on the rapidity of one jet, but the gain depends on the value of the cut on pt of the jets 
and is accompanied by some loss of event rate. As long as hard matrix elements are used to 
generate the pt distributions of the jets, all three methods yield similar event samples. On 
the other hand, significant differences seem to result if one uses a parton shower method to 
generate the if + 2 jet background 3^. In our view, the hard matrix element approach is 
a more faithful representation of the momentum distributions of jets in the relevant WBF 
region of large pt- 



V. ESTIMATES OF THE Hjj NLO CORRECTIONS 

A differential calculation of the Hjj background does not exist at next-to-leading order. 
Prior to undertaking such an effort, we wish to obtain plausible estimates of the sizes of NLO 
effects on both event rates and kinematic distributions. We present two such estimates in 
this section. 

The NLO corrections for QCD production of a Z boson in association with two jets are 



known 



Representative diagrams for this process - seen in Fig. can be compared 
with those for H plus two jets, shown in Fig. El Our first estimate of NLO effects for Hjj 
is based on its similarity with Zjj, but we acknowledge some important differences. In the 
Hjj process the Higgs boson couples only to gluons (via a top quark loop), whereas the 
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FIG. 11: Representative diagrams for leading order Z + 2 jet production. 

Z boson couples only to quarks. This difference means that the processes have a different 
sensitivity to the parton distribution functions. Second, the couplings of the scalar Higgs 
boson to the decay products are also different, and the angular distributions of the decay 
products differ in the two cases. While we do not include a decay branching ratio, our cuts 
require the rapidities of decay products of the produced boson to lie between those of the 
two tagging jets. Finally, the effective Hgg coupling contains a factor of as so that the 
Hjj process is formally proportional to a^, in contrast with for Zjj. Although the Hjj 
process is naively of 0{ag), our use an effective coupling for the Hgg vertex means that 
the Hjj process is effectively of the same order, 0{a1), as the Zjj process. Since we are 
interested only in Higgs boson masses that satisfy_m/^ < 2mf, with transverse momenta 
plf < rrit, the effective coupling approach is valid [2^1. 

We calculate the Zjj cross section using a variable Z-mass, mz = mn- In this way, we 
estimate the NLO corrections for a process involving the QCD production of a heavy vector 
boson and two jets. We examine the distribution of the tagging jet pseudo-rapidities, as 
before. In particular, we are interested in whether the shape of this distribution changes 
significantly at NLO. A NLO effect that would modify the background distribution so that 
it resembles the signal peak, could have a serious effect on the ability to select the genuine 
WBF events. 

We show the lowest order distributions for the Hjj and Zjj processes in Fig. El We do 
not reproduce the plots for more than one value of the pt cut, since this cut does not alter 
the conclusions. One can see that the distributions for the two processes are very similar in 
shape in the two cases, differing in the behavior at low pseudo-rapidities - where the Zjj 
curve is somewhat higher - and toward the tail, where the Hjj distribution dies off slightly 
more slowly. For the purposes of this study, the most prominent difference - that at low 
values of |?7| - has little effect since our cuts require that \rij\ > 1.6 for at least one of the 
tagging jets. The remaining difference in shapes is small, and we conclude that, as a first 
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FIG. 12: Normalized LO tagging jet pseudo-rapidity distributions in Hjj (solid) and Zjj (dashed) 
events. The minimum jet px cut is 40 GeV; rriH = 115 GeV (left) and mn = 200 GeV (right). For 
the Zjj events, the Z mass has been altered to mz = mu- 



PT cut [GeV] 


20 


40 


80 


ruH = 115 


1.00 


0.99 


0.84 


m-H = 200 


1.12 


1.11 


1.02 



TABLE VI: The /C-factors, as defined by Eq. El for Zjj production with our WBF cuts. The Z 
mass has been altered to take on the two relevant Higgs boson mass values, mz = rnn, and the 
scale choice is /u = uih- 



estimate, our use of the Zjj process to approximate Hjj is reasonable. 

For the Zjj process, the lowest order and NLO distributions in the tagging jet pseudo- 
rapidity are shown in Fig. ^1 These are to be compared with the Hjj curves in Fig. El We 
show curves for two different boson masses, as before, and for a variety of pr cuts. The NLO 
corrections do not appear to alter the shape of this distribution significantly. Moreover, with 
the scale choice fi = mn, the corrections are small in magnitude, both in the total cross 
section and over the pseudo-rapidity range of interest. The change in cross section is shown 
quantitatively in Table IVIl where we show the i^-factors for this process having applied all 
the WBF cuts. The corrections vary from ^ 10% for mz = 200 GeV and moderate pt = 20, 
40 GeV to ^ -15% for niz = 115 GeV and high pr = 80 GeV. 

In the definition of the i^-factor in this paper, different parton distribution functions 
(PDF's) are used in the numerator and denominator, LO expressions in the denominator 
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FIG. 13: Dependence of the tagging jet pseudo-rapidities in Zjj events on the minimum jet pT 
cut, for the two cases mz = 115 GeV (left) and mz = 200 GeV (right). In each graph, the LO 
(dashed) and NLO (sohd) curve is shown. The scale choice is fi = mn- 



and NLO expressions in the numerator. Correspondingly, different values of as are used in 
the numerator and denominator: 

^AfLO(CTEQ6M;a™(/i)) 



K 



(5) 



(T^0(CTEQ6Ll;a,^0(/i)) 
That the ii'-factors are close to unity for the Zjj process results from a compensation 
between the change in PDF from LO to NLO and the change in as{n), plus the effects of 
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Pt cut [GeV] 


20 


40 


80 


Signal {niH = 115) 


1374 


789 


166 


Bkg 


1196 


378 


77 


Purity 


0.53 


0.68 


0.68 


Signal {mn = 200) 


928 


545 


121 


Bkg 


598 


199 


47 


Purity 


0.61 


0.73 


0.72 



TABLE VII: Event rates for the Hjj WBF signal(NLO) and Hjj background(estimated NLO), 
for an assumed 1 Po~^ of luminosity. No branching ratio is included for the Higgs boson decay. 



the additional processes at NLO. In as{fi), the net change, after reduction in the value of 
as{Mz) and the altered evolution, tends to decrease the cross section from LO to NLO, by 
a factor 

for both cases (/i = 115, 200 GeV). One must be careful to apply this i^-factor consistently 
only to a lowest order calculation with the same PDF set and treatment of a^. Since we use 
CTEQ6L1 in the background calculation of TablelTTl it is straightforward to incorporate these 
fC-factors in order to obtain the new background estimates in Table I VI II for an assumed 
1 fb^^ of integrated luminosity. For the estimate of NLO corrections to the background 
presented in this section, with /i = mn, the table shows that the purity starts at about 50% 
at rriH = 115 GeV, if the cut on pt is 20 GeV, and grows to about 70% when the pt cut is 
40 GeV or larger. Slightly larger values are obtained at niH = 200 GeV. 

We show in Fig. El the effect of a lower scale choice /i = mn/^ on the Zjj process. In 
contrast to Fig. E[ the NLO corrections are now substantial, negative, and not constant 
as the pseudo-rapidity changes. For this lower scale choice, we now apply the /^-factors 
point-by-point to the lowest order Hjj background distribution in order to estimate the 
new NLO result. The result of this exercise is shown in Fig. along with the estimated 
NLO result for ^ = mn- This figure shows that our estimate for the NLO Hjj background 
cross section is affected very little if we choose a smaller scale such as /i = rriH/^. Reducing 
the scale, we find that the LO background is increased substantially. However, the i^'-factor 
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2 4 2 4 



FIG. 14: Tagging jet pseudo-rapidity distributions in Zjj events, calculated with a smaller renor- 
malization and factorization scale, ji = In each graph, the LO (dashed) and NLO (solid) 

curve is shown. 




FIG. 15: Estimated tagging jet pseudo-rapidity distributions in Hjj events at NLO, calculated 
with a smaller renormalization and factorization scale, // = mj|//2. Also shown (lower curve) is the 
estimated NLO Hjj result with /i = mn- 

decreases in such a way as to restore the size of background at NLO. 



A. Background estimate with a mass cut 

In a similar spirit to the above study, here we present the estimate of the NLO corrections 
to the Hjj background when the definition of the WBF event sample involves only an 
invariant mass cut on the two tagging jets, Mjj > 800 GeV. In Fig. UHl we show the lowest 
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Pt cut [GeV] 


20 


40 


80 


tuh = 115 


0.89 


0.80 


0.78 


mu = 200 


0.93 


0.93 


0.89 



TABLE VIII: The /C-factors, as defined by Eq. for Zjj production with an alternative definition 
of the WBF cuts, Mjj > 800 GeV. The Z mass has been altered to take on the two relevant Higgs 
boson mass values, mz = rnn- 



order distribution of the invariant mass of the two tagging jets, for the processes Hjj and 
Zjj and a jet pt cut of 40 GeV. Each curve is normaUzed by its own integrated cross section. 
Although the two distributions differ in shape considerably over such a wide range of Mjj, 
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FIG. 16: Normalized LO tagging jet invariant mass distributions in Hjj (solid) and Zjj (dashed) 
events. The minimum jet pT cut is 40 GeV; niH = 115 GeV (left) and niH = 200 GeV (right). For 
the Zjj events, the Z mass has been altered to mz = ran- 



the two curves show similar behavior above Mjj > 800 GeV. We conclude again that the 
Zjj process will yield a reasonable estimate of the Hjj process in the region of phase space 
of interest for WBF studies. 

The effects of the NLO corrections for the Zjj process are shown in Fig. El They are 
small over the entire mass range, and they do not change the shape of the distribution. 
The net effect is summarized in Table IVllll where we show the i^'-factors for each cut 
and Higgs mass for the alternative definition of the WBF sample, Mjj > 800 GeV. The 
corrections are universally negative, ranging from ^ —10% for mz = 200 GeV to ~ —20% 
for mz = 115 GeV and the higher pt cuts of 40,80 GeV. 
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Mjj [GeV] Mjj [GeV] 



FIG. 17: Tagging jet invariant mass distribution in Zjj events with a jet pT cut of 40 GeV, for 
the two cases mz = 115 GeV (left) and mz = 200 GeV (right). The LO (dashed) and NLO (solid) 
curves are shown, for a scale choice of /x = mn- 



PT cut [GeV] 


20 


40 


80 


Signal {rriH = 115) 


808 


561 


158 


Bkg 


271 


146 


64 


Purity 


0.75 


0.79 


0.71 


Signal {mn = 200) 


617 


428 


121 


Bkg 


146 


88 


38 


Purity 


0.81 


0.83 


0.76 



TABLE IX: Event rates for the Hjj WBF signal(NLO) and Hjj background (estimated NLO), 
without our usual WBF definition and instead with Mjj > 800 GeV. We have assumed 1 fb^^ of 
integrated luminosity. 

Using these i^'-factors, we can update the earher leading order results of Table M Since 
the i^'-factors are close to unity, the effect of the NLO corrections is small, as illustrated in 
Table HXl With NLO effects included, we conclude as at LO, that the cut on the jet-pair 
invariant mass improves the signal purity significantly with respect to our simple rapidity 
prescription only for the smallest of the cuts on px {c.f. results Table EH}, but at a cost in 
the signal rate at all pt- 
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B. Second Estimate of NLO Corrections to the Background 



A different and larger estimate of tlie NLO corrections to the Hjj background may be 
obtained if we begin witli results for inclusive Higgs boson production processes. Corrections 



to the fully inclusive cross section are known both to NNLO 2^-j28| and with all-orders 
soft-gluon resummation included [s^-js^. In addition, the NLO corrections to the H + j 
process are published js^. By starting from the inclusive process or the H + j process, we 
begin with a process that is dominated by gg scattering, unlike the Zjj case. 

With the same definition of the i^'-factor used in this paper, and the same choice of 
renormalization/factorization scale, the result for the fully inclusive Higgs boson cross section 
is 

K{H + X) = 1.7-1.8. 

The increase at NLO is slightly larger at rriH = 200 GeV than at rriH = 115 GeV. For the 
semi-inclusive H + j process, where an additional jet of pt > 20-80 GeV is specified in the 
final state, the NLO corrections correspond to, 

K{H + j+X) ^ 1.3 - 1.5. 

In this case the PDF sets used were similar but not identical (CTEQ4), and this result is 
quoted for mn = 120 GeV. 

We now try to apply these K factors, obtained in situations in which either no jet or at 
most one jet is required in the final state, to the case of interest here where two jets are 
required, each with a minimum value of pt, and with the other WBF cuts. In the case 
of Z + jets, the i^-factors are 1.14, 1.16, and 0.90 for 0, 1, and 2 jets, respectively {2^, 
with the same PDF's and definition of i^-factors used here. In other words, the i^'-factor 
is approximately independent of the number of specified jets. In Table IVIf ii'-factors are 
shown for the Zjj process with WBF cuts applied. They are close to unity. In Ref . 2^ , the 
i^-factor is obtained for the Z + 2 jet process with a jet cut of pt > 20 GeV, but without 
the WBF cuts. It is also close to unity, K = 0.9. This value is similar enough to the 
numbers shown in Table IVTl to justify the assumption that the overall inclusive ii"-factor is 
not very different from the one that should be applied for WBF events. Adopting this line 
of argument, we suggest that a conservative estimate of the i^'-factor for Hjj production in 
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Pt cut [GeV] 


20 


40 


80 


Signal {niH = 115) 


1374 


789 


166 


Bkg 


1914 


611 


123 


Purity 


0.42 


0.56 


0.57 


Signal {mn = 200) 


928 


545 


121 


Bkg 


854 


286 


74 


Purity 


0.52 


0.67 


0.62 



TABLE X: Event rates for the Hjj WBF signal(NLO) and Hjj background (estimated NLO, 
according to the procedure in Section IV B|) . for an assumed 1 fb^^ of integrated luminosity. 

the WBF region is a i^'-factor from the more inclusive processes above, 
K(„ + + A-) « '<(H + X) + K(H + ,^X) ^ 

The new estimates of NLO corrections to the background are presented in Table |X| The 
signal purity starts at about 40% when mn = 115 GeV, if the cut on px is 20 GeV, and 
grows to about 60% when the pt cut is 40 GeV or larger. Slightly larger values are obtained 
at tuh = 200 GeV. 

With the larger background estimate of this subsection, the signal purity is lower than 
before (c./. Table 1X1 vs Table rvTT]) . If the NLO enhancement of the background is as great 
as a factor 1.6, the value of P ~ 0.4 presented in Table IX| suggests that investigations with 
a low Pt cut will be ineffective. According to the results in Fig. [3 the uncertainty in 6g/g is 
quite large when P < 0.5. With the cut pt > 40 GeV, the corresponding purities of P ~ 0.6 
may still permit determinations of 6g/g ~ 10%. A fully differential NLO calculation of 
the background process for + 2 jets is definitely needed to establish both the size of the 
background and the theoretical uncertainty 6B/B. 

VI. SUMMARY AND CONCLUSIONS 

In this paper, we investigate the weak boson fusion process for production of the neutral 
Higgs boson H at the LHC. We are interested in estimating the accuracy with which the 
Higgs boson coupling to weak bosons may be determined from data. An important and, 
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in fact, controlling aspect is the extent to which events produced by the WBF subprocess 
may be separated from events in which a Higgs boson is produced by other mechanisms. A 
hallmark of the WBF subprocess is that the Higgs boson is accompanied in the final state 
by two jets that carry large transverse momentum and relatively large rapidity. However, 
purely strong interactions subprocesses also produce Higgs bosons accompanied by two jets. 
To extract the couplings reliably, a good understanding is required of both the production 
and the background processes. We use hard QCD matrix elements in order to represent 
the signal and the + 2 jet background reliably. We provide an independent calculation 
that verifies the fully differential next-to-leading order QCD corrections to the WBF signal 
process of Ref. ^|, and we examine in more detail the effects of the WBF selection cuts 
on these NLO QCD corrections. We use leading order perturbative QCD expressions for 
the background if + 2 jet matrix elements since NLO results are not yet available in fully 
differential form. We also provide two estimates of the NLO enhancement of the QCD 
if -|- 2 jet background process. Our calculations are fully differential at the partonic level. 

Among our goals in this study are to evaluate the effectiveness of different prescriptions 
for defining the WBF sample and to estimate the expected WBF signal purity P, by which 
we mean the fraction of real Higgs boson events produced by weak boson fusion. 

Various prescriptions are used in the literature to define the WBF sample, cuts that en- 
hance the WBF fraction of the cross section by exploiting the special transverse momentum 
and rapidity characteristics of WBF events. Our investigations lead us to propose a new, 
somewhat simpler definition in terms of a cut on the rapidity of one of the final state jets, as 
defined by Eq. (^. We compare the effects on both event rates and signal-to-background ra- 
tio of our proposed method for defining WBF events with two other popular methods found 
in the literature: a selection on the difference in rapidities between two tagging jets in the 
final state, and a selection on the invariant mass of a pair of tagging jets. In a low-luminosity 
environment where one may be tempted to use a relatively low cut on transverse momentum 
to select the trigger jets, the conventional alternatives provide better signal purity but at a 
cost in signal rate. Once the cut is raised, our definition does essentially as well in signal 
purity while preserving more of the signal. 

We find that purities of 60% to 70% can be expected if a selection of pt > 40 GeV is 
made on the tagging jets. We derive an expression for the expected uncertainty on the effec- 
tive Higgs-boson-to-weak-boson coupling strength g in terms of P, the expected statistical 
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accuracy of the LHC experiments, and the theoretical uncertainties on the signal and the 
background processes. We estimate that an accuracy of 5g/g ~ 10% may be achievable 
after ~ 200 fb~^ of integrated luminosity is accumulated at the LHC. On a cautionary note, 
however, we recall that our WBF signal purity and our uncertainties are obtained in a very 
well controlled situation in which there is an identified Higgs boson in a sample oi H + 2 jet 
events produced by both the WBF mechanism and the QCD background processes. In an 
experimental context, there will be additional sources of background from final states that 
mimic a Higgs boson. The effects of these additional backgrounds presumably only increase 
the expected uncertainties on the couplings. 

The theoretical uncertainties on the signal S", and on both the size and uncertainty of 
the background B dominate the uncertainty in g. Current estimates of 5S/S are in the 
5% range, and, since differential NLO calculations exist, this uncertainty is controlled by 
uncertainties in the parton densities and by the residual renormalization and factorization 
scale dependence. In order to reduce the estimated uncertainty in g, the next major step 
would appear to be a fully differential NLO calculation of the background applicable in the 
region of interest for WBF investigations. 
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